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Abstract: It is shown that large Nc QCD must have a Hagedorn spectrum (i.e., a spec- 
trum of hadron which grows exponentiahy with the hadrons' mass) provided that certain 
technical assumptions concerning the applicability of perturbation theory to a certain class 
of correlation functions apply. The basic argument exploits the interplay of confinement 
and asymptotic freedom. 
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1. Introduction 

More than four decades ago, Hagedorn ^ proposed that the number of hadrons with mass 
less than m grew exponentiahy with m. Neglecting hadron-hadron interactions this implies 
that the energy density diverges at a finite temperature, yielding a maximum temperature 
for hadronic matter, Th- A Hagedorn spectrum was found to arise automaticahy in string 
theories ^ (which were originahy formulated as a theory of the strong interaction). Thus, 
the idea of the Hagedorn spectrum fits naturahy with our understanding of QCD: QCD is 
confining and as such one might expect it to have stringy dynamics for high-lying states. 
Moreover, at high temperatures QCD is in a quark-gluon phase rather than a hadronic 
phase. 

Despite empirical evidence that the number of hadrons up to ~ 2 GeV does grow 
rapidly with the mass in a manner qualitatively consistent with an exponential growthQ, 
the picture described above is clearly imprecise. Hadron masses are not strictly well defined: 
hadrons have widths due to strong decays. Moreover, unlike in pure Yang-Mills theory, 
there is no order parameter for deconfinement in QCD; as such, the distinction between a 
hadronic phase and a deconfined phase is not sharp. However, the notion of a Hagedorn 
spectrum becomes precise if applied to the large Nc limit of QCDQ, ^] rather than QCD 
itself. In the large Nc limit meson and glueball widths go to zero making the hadron 
masses well defined[Q, Moreover, at large Nc, QCD has a first order phase transition to 
a deconfining phase [^. 

The conventional wisdom is that large Nc QCD has a Hagedorn spectrum: the number 
of hadrons with mass less than m (denoted N{m)) in large Nc QCD satisfies the condition 
that for sufficiently large m 

N{m) > e™/^ (1.1) 
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Th is defined as the maximum value of T for which inequahty ( |1.1[) holds. It is known 
that large Nc QCD with adjoint fermions in 1+1 dimensions has a Hagedorn spectrum^. 
While there is no real reason to doubt a Hagedorn spectrum also holds for large Nc QCD 
in 3+1 dimensions, to date there has been remarkably little direct evidence for this. There 
has been a study of the large Nc glueball spectrum based on a numerical treatment of a 



transverse lattice in a light cone formalism] 11 1. The results were consistent with a Hagedorn 



spectrum-but not definitive. Indirect evidence for a Hagedorn spectrum at large Nc can 
also be obtained from lattice studies of thermodynamic properties of hadron matter in the 
metastable above the deconfining transition [p!o|]. 

The purpose of the present paper is to sketch a first principles argument that a Hage- 
dorn spectrum must arise in QCD (i.e., Yang-Mills fields plus quarks in the fundamental 
representation) at large Nc- The approach relies only on generally accepted properties 
of QCD such as asymptotic freedom and confinement (in its basic sense that all physical 
states are color singlets) and on plausible assumptions about the applicability of pertur- 
bation theory to describe correlation functions at short distances. It explicitly assumes 
neither that confinement is manifest through an unbroken center symmetry nor that the 
dynamics of the hadron is stringy. The approach is based on the fact that the number of 
independent local operators with fixed quantum numbers and a single color trace grows 
exponentially with the mass dimension of the operator. This approach is similar in spirit 
to the demonstration by Kogan and Zhitnitsky that large Nc QCD with adjoint fermions in 
1+1 dimensions has a Hagedorn spectrum It also has elements which are reminiscent of 
refs. 1^, P]; these works deal with a rather different problem, namely, the thermodynamics 
of systems confined to a small sphere rather than spectroscopy. 

The approach relies on a confrontation between asymptotic freedom with confinement 
and depends on working in a regime in which perturbative corrections to the leading free- 
field values for certain correlation functions is small. At first glance it may seem to be 
impossible to learn anything about the hadronic spectrum directly from the perturbative 
regime: the perturbative regime by construction only describes the smooth part of spectral 
functions. However, this is misleading. While one cannot learn anything directly about 
individual hadrons, one can learn certain global features about the spectrum, at least at 
large Nc- Recall that in the formal large Nc limit mesons are narrow. High-lying mesons 
do not "melt" into the continuum due to large widths associated with phase space for 
decay. Thus the spectral functions for correlators with meson quantum numbers are given 
entirely by mesons poles up to arbitrarily high mass. The perturbative regime is for the 
deep space-like region of momentum space far from these poles; there is no ability to resolve 
information about individual hadrons in this regime. However, dispersion relations relate 
the space-like correlators to integrals over the spectral functions which are given by the 
meson poles. Thus, while the perturbative regime cannot tell us about individual hadrons, 
it does contain some global information about the high-lying spectrum at large Nc- As will 
be shown in this paper, the existence of a Hagedorn spectrum is global information about 
the spectrum which is extractable, provided one accepts certain technical assumptions 
about the applicability of perturbation theory. 

Before turning to the Hagedorn spectrum, it is useful to recall a classic example where 
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the nature of the asymptoticahy free regime for correlators is used to infer basic properties 
of the high-lying spectrum: the demonstration that large Nc QCD has an infinite number 
of mesons with any given quantum numbers {eg., the scalar-isoscalars) Q. Witten's orig- 
inal version of this argument was based on momentum space correlators. Here it will be 
slightly recast in terms of correlators in (Euclidean) position space. This is useful since the 
argument for the existence of a Hagedorn spectrum is based on position space correlators. 
Consider the correlator of the scalar current J = qq at two different space-times points. 



By the standard Kalleh-Lehmann spectral representation! 12 1 it can always be written as 
n(r) = (J(f)J(O)) = / dsp{s)A{r;s) 



with A(r; ^) ~ J 



(27r)4 g2 + s 27r2r dr 



where r = + and p{s) is the spectral function. Large planarity implies that 
all physical states created by the current is composed of quarks and gluons in a single 
indivisible color singlet combination (plus ^jNf. corrections) [Q, |5|. If one imposes on this 
structure confinement (in the very basic sense that all physical states are color singlets), 
then one deduces that at large the currents can only makes single meson states [Q, |5|. 
The spectral function at large is thus given by 

p{.s) = Y,\c,fKs-rr?^ (1.3) 
i 

where vrtj is the mass of the j*'^ meson. 

The scalar propagator, A, has the property that 

^('^> Tr4-T as r ^ (1.4) 

for any s; this asymptote is approached when r\fs <^ 1. Now suppose there were only 
a finite number of scalar-isoscalar mesons, then there must be a meson with a maximum 



mass. In that case, Eqs. (1.2) and (|5.2|) together with Eq. (1.4) implies that 



Jmax 



n(r) ^ |cilV2vrV2 (1.5) 
i=i 

for r <C "^j^^^- However, asymptotic freedom implies that when r ^ (where A is 
the QCD scale), H ~ r~^. This is inconsistent with the scaling of Eq. ( [1.51) . Thus we 
conclude that the assumption that there are only a finite number of mesons used to derive 
Eq. ( |1.5| ) must be false: there is an infinite number of mesons at large Nc- 

A couple of comments about this example are in order, first, the example shows that 
knowledge of correlation functions in the perturbative region is sufficient to extract some 
qualitative information about the high- lying spectrum of mesons for large Nc QCD. The 
purpose of this paper is to show that the existence of a Hagedorn is similarly a qualitative 
property of the high-lying spectrum which is accessible from information in the perturbative 
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regime, second, it should be apparent that the procedure used in this example depended 
critically on taking Nc ^ oo prior to the small r limit. It was only by taking A'^c — ^ oo at 
the outset that one could justify the use of Eq. ( ^.2| ) for small r. In what follows, it will 
always be assumed that the large Nc limit is taken prior to any others. 

The approach to establishing Eq. (^]^) is to generalize J from a single operator to 
large sets of linearly independent local operators all with the quantum number of interest. 
Rather than having a single correlator as in the preceding example, one studies a matrix 
of correlators based on a set of local operators. 

The critical fact underlying this approach is the fact that the number of distinct oper- 
ators with fixed quantum numbers grows exponentially in the dimension of the operator. 
This is easily shown. Ultimately, this will translate into the condition that the number 
of distinct hadrons at large Nc must grow exponentially with the mass provided that the 
regime of validity of perturbation theory is as expected. Note that while there is a connec- 
tion in this derivation between the exponential growth of the number of distinct operators 
with dimension and distinct hadrons with mass, the connection is somewhat subtle. The 
derivation does not assume that each operator couples to distinct hadronic states; and 
indeed they do not. Instead the derivation exploits an inequality relating the lowest mass 
hadrons states with fixed quantum numbers to properties of the trace of the logarithm of 
the correlator matrix. 

The strategy for deriving a Hagedorn spectrum has four basic components. The first 
is based on the fact that as the spectral functions for the correlator matrix are saturated 
by meson poles at large Nc (a fact also used to prove that there were an infinite numer 
of mesons at large Nc). From this, plus general properties of the scalar propagator, one 
can derive an inequality bounding the average of the lowest k meson masses (where k is 
the number of operators) from above by the derivative (with respect to r) of the average 
diagonal element of the logarithm of the matrix of correlators. 

The second is a generalization from a single set of operators to an infinite sequence 
of sets of operators: Si, S2, ■ ■ .... The key thing is that one can show that the 
inequality mentioned above implies a Hagedorn spectrum provided the dimension of the 
set of operators grows exponentially in n (the index specifying the term in the sequence) 
{i.e., log(fe) ~ n), while the derivative (with respect to r) of the average diagonal element 
of the logarithm of the correlator matrix for the set of operator asymptotes (at small r) 
to a constant times —n/r with a correction whose relative size is bounded in a particular 
way. 

The third ingredient is an explicit construction of such a sequence of sets of operators 
which one expects to meet the above criteria. Here the focus will be on operators with 
the quantum numbers of scalar and pseudoscalar mesons. The construction begins with 
local gluon bilinear operators coupled to a Lorentz scalar but in the color adjoint. There 
are two such operators; one with positive parity, the other with negative parity. The rfi^ 
element of the sequence is the set of all single-color-trace operators composed of an anti- 
quark operator followed by n of these bilinears followed by a quark operator. This set 
clearly grows exponentially as there are 2" elements. Moreover, at asymptotically small r 
and large n asymptotic freedom requires that derivative of the average diagonal element of 
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the logarithm of correlator matrix asymptotes to —8n/r, satisfying two of the conditions 
needed to demonstrate a Hagedorn spectrum. 

The final step is to show that the correction to 8n/r has a relative size that is bounded 
appropriately. At this point it is necessary to rely on perturbation theory. It can be shown 
that perturbative corrections at any fixed order satisfy the requirement. The critical point 
is to show that the correction at any given order in perturbation theory has a contribution 
independent of n, plus 1/n corrections but no terms which grow as n. Thus, to the 
extent that perturbation theory is valid — as is generally expected for correlators at short 
distance — a Hagedorn spectrum must emerge. 

It is clear from this brief description, that argument does not constitute a rigorous 
theorem. The reliance on perturbation theory is problematic. In the first place it is 
approximate and thus difficult to use setting a precise bound, second, it is an asymptotic 
series. Due to the asymptotic nature of the series it is very difficult to see how to tighten 
this argument into a theorem. The usTial assumption in QCD is that perturbation theory 
announces its own demise. That is to say that if one studies correlation functions beginning 
in the asymptotically free regime of short distances or high momenta, and then pushes to 
longer distances or lower momenta perturbative corrections (calculated at some order in 
perturbation theory) grow. To the extent that these corrections remain a small fraction of 
the total, the usual expectation is that perturbation theory should be valid. The derivation 
here depends on this expectation holding for the matrix of correlation functions in question. 
While it should be stressed that this is an assumption and has not been proven rigorously, 
it should also be stressed that this assumption is quite standard. To the extent that this 
expectation is correct, a "physicist's proof" of a Hagedorn spectrum might be said to exist. 

The fact that argument is less than completely rigorous is hardly surprising — virtually 
nothing about QCD is known with full mathematical rigor including such basic features as 
asymptotic freedom which is also demonstrated perturbatively. However, the argument is 
of value for several reasons, despite its lack of rigor. First, there is considerable experience 
that perturbation theory does accurately describe short-distance correlators; to the extent 
that this is true the Hagedorn spectrum has been established. Second, the argument can 
be easily generalized to various types of QCD-like theories with different numbers and 
types of fermions, different dimensions of space-time and different gauge groups. The 
approach shows which of these have Hagedorn spectra and why. Third, the argument 
does not explicitly assume that QCD approaches some type of a string theory for highly 
excited states. To the extent that one takes a Hagedorn spectrum as a signature for stringy 
dynamics this argument gives insight into the emergence of stringy dynamics while only 
using generic properties of correlators. Finally, the argument may provide some insight into 
the nature of confinement. The argument only uses confinement in the sense that physical 
states arc color singlets; but it does not explicit rely on confinement being manifest through 
an unbroken center-symmetry. 

This paper is organized as follows: the next four sections will discuss each of the four 
basic components of the argument. Following this will be a discussion on how the argument 
can be generalized to various QCD-like theories at large Nc and some closing remarks. 
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2. A useful Lemma 



The first step in estabhshing Eq. (LI) is to generahze J from a single current to large 
sets of linearly independent local gauge invariant operators constructed from quarks and 
gluons. These operators will all have fixed quantum numbers of interest and all involve only 
a single color trace. This last condition ensures that at large Nc the operator, when acting 
on the vacuum only creates single hadrons. For technical simplicity here, these operators 
will be restricted to those with the quantum numbers of scalar-isoscalar mesons (where 
scalar refers to spin but not parity; pseudoscalars are included). It is straightforward to 
generalize the argument to show exponential spectra for mesons of other quantum numbers 
and to glueballs, but the argument is particularly straightforward for scalar mesons. 

In what follows, an arbitrary constant with dimensions of mass, A is introduced and 
inserted in appropriate places so that correlation functions are dimensionless. The precise 
value of A is irrelevant as A cancels out in all final results. 

Begin with a set of k linearly independent currents, single color trace, scalar-isoscalar 
currents S = {Ji, J2, J3 ' ' ' Jk}- Associated with 5 is a dimensionless kxk mixed correlator 
matrix ^^{r). For any operators Ja,Jb^S the matrix elements are given by 

n-'-W ^ = / dspab{s)/\{r;s) with /,,,(s) = c;,.c,,,<^(s - m|); (2.1) 

j\ is the Hermitian adjoint of J a , dc is the mass dimensions of operator and is 
the (dimensionless) amplitude for the a**^ current to create the j^^ scalar-isoscalar meson: 
Ca,j = A~'^'^(j| Ja|vac) (where \j) is the hadronic state associated with the j*'^ scalar-isoscalar 
meson). The form of the spectral function follows from the usual large Nc requirements. 

Given a set of such operators, 5, an important lemma can be established for meson 
masses in the large Nc limit. The lemma is that for any r > 0, 

^ tr(log(1f^(r))) I 



mi is the lightest scalar-isoscalar meson, m2 the second lightest, etc., and 1 15| | = dim (^^^[r] 
k. Note that this lemma bounds the average of the lowest k mesons by minus the derivative 
of the average diagonal matrix element of the logarithm of the correlator matrix. 

This lemma can be established from general properties of linear algebra, Eq. ( p.l] ) 
and along with a few basic properties of the scalar propagator function as discussed in 
Appendix |^. 



3. A sequence of sets of operators 

Consider a sequence of sets of scalar-isoscalar linearly independent gauge invariant local 
operators with quantum numbers of a scalar-isoscalar meson and composed of single color 
trace operators, 5i, ^2, ^3, . . . , 5„, . . . . Suppose that such a sequence exists with the 
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property that one can always find positive constants rg, A, p, and uq such that for all 
r < ro and all n > uq the following two conditions are satisfied: 



\Sn\\ > A" ; b) n^>MsAr) ■ (3.1) 
r 



Given such a sequence one can establish the Hagedorn spectrum of Eq. (1.1). 

To see this consider A'o(m), the number of scalar-isoscalar mesons with mass less than 
m. By construction N{m), the total number of hadrons with mass less then m, is greater 
then or equal to NQ(m). Define 

where rrij is the mass of the jth scalar meson; i.e., (m)o is the average mass for all scalar 
mesons with mass less than rriQ. Start by assuming that No{m) grows with m at least 
linearly. In that case (m)o > mQ/2. Choose a value of uiq such that NQ{mo) = 
Using the lemma in Eq. (|2.2[) and condition b) in Eq. (|3.1|) allows one to deduce that 



Thus 2np/rQ > niQ. Acting on both sides of this with the monotonic function Nq yields 

No{2np/ro) > No{mo) = \\Sn\\ > A"" (3.4) 



where the last inequality follows from condition a) of Eq. ( |3.1j ). Finally, denoting 2np/rQ 
as m yields 

N{m) > No{m) > (A'^o/^p)'" = g'nCro iog(A)/2p) ^3^5^ 
This is of the form of a Hagedorn spectrum with 

Th < (3.6) 
ro log(A) 

This derivation depended on the assumption that No^m) scaled at least linearly. This 
assumption is, of course, justified a posteriori by the form of Eq. (|3.5|) . While the preceding 
argument is sufficient to obtain a Hagedorn spectrum, it is worth noting that the bound 
on Tff in Eq. ( p. 6]) can be sharpened. Note that the factor of 2 came from the assumption 
that NQ{m) grew at least linearly. If instead the assumption had been that Nq grew at 
least as fast as a given power law: N{m) > const x A^" (which is also justified a posteriori) 
then the condition on Th becomes Th < a roTog^( A) Taking a to be arbitrarily 

large yields 

Th < . (3.7) 

rolog(A) 
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4. Constructing operators 



The next step is to construct explicitly a sequence of sets of operators which satisfy the 



conditions in Eq. (3J) and thereby establish a Hagedorn spectrum. 

Two useful ingredients in this construction are the scalar and pseudoscalar operators 
composed of bilinears in the gluon fields. 



0+^^F,,F^'^ O^^^F,,F^^. (4.1) 



These operators are not traced over color; and at large they become pure color adjoint 
operators up to 1/Nc corrections. From these one can construct operators of the following 
form 

Jl,,l^,...^l^ = qOi,Oi^---Oi^q where Zi,--- =± (4.2) 

eg., J+_ = qO+O-q = qFf,^F'"'Fc,f}F°''^q/N^. Define 5„ to be the set of ah operators of 
this form of (engineering) dimension 4n + 3, i.e., containing n Oi operators. Thus, 



5i = {J+,J_} 



'qF\^^^F^ qF\^J^\ 



= {J,,, J,_, J__} = \ ^^^-^Z"^^"'' . (4-3) 



iV4 

^3 = {J+++, J++-, J+ ) J-+-, • • • } 



By construction = 2" so that the sequence 5i, 52, ^3, . . . , 5„, . . . obviously satisfies 



condition a) of Eq. (3A); demonstrating that condition b) holds is equivalent to establishing 



a Hagedorn spectrum. 



5. Establishing the conditions needed for a Hagedorn spectrum 

To establish condition b) of Eq. ( |3.lD for this sequence of sets of propagators — and thereby 
establish the existence of a Hagedorn spectrum — one begins with an analysis of the short- 
distance behavior of ^?(r). As r — )• 0, asymptotic freedom implies that correlators take 
their free-field values (as represented for example by diagram a. of Fig. |l|). Since the two- 
gluon state created by 0+ associated with some particular color indices are free fields they 
cannot be annihilated by 0-. Moreover, the chance that the same color index appears 
more than once among the 0± operators vanishes at large Nc- Thus at large Nc and r — t- 
the correlator will vanish unless all of the 0± operators in each current occurs in the same 
order: ^ is diagonal. 

Moreover when r is much smaller than Aq^j-,, simple dimensional analysis dictates the 
form of the diagonal matrix elements of tf ; they must all be of the form const r-(««+^) 
since is of mass dimension 8n -|- 6 and r is the only dimensional parameter left at those 
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d. 



e. 



f. 
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Figure 1: Typical diagrams which contribute to diagonal correlators associated with the currents 
in iS". The solid lines correspond to quark propagators. The dotted lines (dashed lines) represent 
the free propagation of two gluons coupled to a scalar (pseudo-scalar) as created by 0+ {0-). 
The black circles correspond to interactions due to one or more gluon exchanges. The smaller 
black circles represent gluon exchanges within one pair of gluons while the larger circles represent 
gluon exchanges between two pairs. The ellipsis indicates pairs of gluons which are not explicitly 
represented. The large grey ovals correspond to operators in 5" . In diagrams a. through f . the initial 
operator and final operators are the same; these diagrams contribute to diagonal matrix element in 
the matrix of correlators. In diagram g. the initial and final operators differ and corresponds to a 
contribution to a mixed correlator, i.e. an off-diagonal matrix element. 

scales. Prom this it is easy to see that Ms,^{r) — )■ 8n -|- 6 as r — >■ 0. Given this asymptotic 
form, it is useful to parameterize M^^ (r) in the follow way 



where i?(r, n) is the size of the correction due to interactions. The condition for a Hagedorn 
spectrum amounts to the condition that there exists a number p such that 




8 + ^ 




r 




whenever n> uq and r < ro ; 




r 



in that case condition b) of Eq. (|3.lD is satisfied with p = ^ + Q/uq + p. 

To study the behavior of i?(r,n), it is useful to parameterize If as the product of its 



free field value times a correction due to interactions: 

(5.3) 



log (tf (r)) = log (tf frce(r)) - ^-Y^'^'ir 
From the definition of Mg^ (r) and R one sees that 



Ideally one should obtain a rigorous bound on R{r,n). It is not immediately appar- 
ent how to do this. However, given that the objects under study are correlation func- 
tions, it is natural to make the standard assumptions used in the study of QCD corre- 
lation functions — namely, that at short distances correlators are accurately described by 
renormalization-group-improved perturbation theory. Under this assumption, it will be 
shown that condition b) holds when working up to any given order in perturbation theory; 
it will then be argued that non-perturbative effects are not expected to alter this conclu- 
sion. Before proceeding with a perturbative analysis it is worth noting that this result 
is nontrivial. It is not sufficient to argue that at short distances perturbative corrections 
generically become small since becomes small. The key issue is that generically, the 
coefficients in a perturbative expansion (i.e., an expansion in ag) for R can be expected to 
depend on n. If any of these coefficients grow with n faster than linearly then condition 



b) of Eq. (3.1) is not satisfied for the perturbative estimate of Ms„- Moreover, the typical 
behavior of combinatoric factors in perturbative theory might lead one to expect that co- 
efficients would diverge with n. However by taking the large Nc limit at the outset, these 
combinatoric factors are greatly suppressed, yielding a linear dependence on n. 

Suppose the correlators in the set of Sn (specified above) are computed perturbatively. 
The expansion at large Nc is in the 't Hooft coupling N^ag- Suppose one works to some 
fixed order, (A^c^s)', where Us is taken to be evaluated at the scale /x^ = r~^. It can then 
be shown for n>l + l that 

R{r, no) = drinY, (A^c".)' 9i + Yl (^-«^)' (^-5) 

\ i=l i=l / 

where gi and hi are numerical coefficients. Note that the only r dependence on the right- 



hand side is through the scale in a^. The two essential features of Eq. (5.5) are that the 
coefficients gi and hi are universal — their values do not depend on n — and that at large n, 
the dominant contributions scale with n but not faster. 

The derivation of Eq. ( p.5p is somewhat involved. However, the origins of its key 
features are relatively easy to understand. One critical ingredient is the use of correlators 
in position space. This ensures that contributions of non-interacting parts of a correlator 
simply factorize: the contribution to a given diagram of a particular type of cluster involving 
the interaction of some number of 0± currents is independent of n. This gives rise to the 
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universality of the coefficients gi,hi. A second key ingredient is the logarithmic structure 
of M5„. 

The final critical ingredient in obtaining Eq. ( |5.5| ) is the large Nc limit. This sup- 
presses non-planar diagrams and thus implies that clusters of interacting gluons only in- 
volve neighboring 0± currents. This limits the combinatoric growth of the amplitudes 
with n — ultimately restricting it no faster than n. To illustrate this, consider one gluon 
exchange contributions to some diagonal matrix element of ^ . The Feynman diagram 
contains 2 quark lines and 2-n gluons which enter and leave the operators. Thus, at the 
one-gluon exchange level there are (2n + 2)^ distinct one-gluon exchange contributions be- 
tween these lines. This grows faster than n. However, most of these gluon exchanges are 
non-planar and thus suppressed at large Nc- There are only 2n + 1 planar contributions 
and this growth is linear in n — as is required. 

The logarithmic structure of Mg^ plays an essential role in restricting the scaling of 
R{r, n) with n to no faster than linearly. To see why, consider diagrams b. and c. of Fig. |l]. 

Diagram b. of Fig. || represents the contribution of all gluon exchanges up to fixed 
order j in perturbation theory contained within one the positive parity pair of gluons. In 
effect it is a self-energy contribution for the pair as a whole — for one pair in one diagonal 
matrix element in ^ (r). Denote this contribution to the diagonal matrix element as S4.. 
Note that S+(r) can be calculated perturbatively and expressed as a series in NcOtg. Next 
note that because the correlator is point-to-point, contributions factorize. Thus one has the 
same self-energy contribution to this matrix element from any positive parity pair; the total 
contribution to the matrix element from all diagrams with a single self-energy of a positive 
parity pair is n+Sj-, where n_i_ is the number of positive parity pairs in the operator J„. 
The LLbutJn oVall such L^.,., to t.- (S (^) . thu' („I)L (where the averaging 
is over states in the class) yielding a total contribution of nS_|_/2; the contribution to 
R{r,n) from the j = 1 term in Eq. ( ^.4| ) from these diagrams is ndrTi^/2. As expected 
this contribution to R{r, n) grows linearly with n. 

However, there are also contributions which scale as n^. First of all there is the 
j = 2 term in Eq. (^.4|); this depends on tr ^^"^{rfj . The contributions to this term from 
diagram b. acting twice is —{{n? + n)/8)5r.S^. There is also a contribution from diagrams 
containing two distinct self-energies such as seen in diagram c . The j = 1 contribution for 
these is {{n? — n)/8)9rS^. Note that this contribution exactly cancels the previous one 
yielding a total contribution to R{r, n) proportional to n. The reason for this cancelation 
is clear: the combinatoric factors are identical. 

This behavior is generic. Cancelations between separate interactions within diagrams 
and higher powers in ^^(r) always occur yielding a linear growth in n and not faster. It 
occurs for interactions involving clusters of adjacent currents. An example is in diagram 
d. of Fig. |l|, which represents gluon exchanges which couple one negative parity pair of 

gluons with a positive parity pair. Denoting the amplitude for this as S |_ is easy to show 

that the contribution from all diagrams containing one such pair to the j = 1 contribution 

to R{r,n) is (n — 1))S ^/4. The contribution of order from j = 2 from diagrams of 

type b. and d. cancel against the order contributions from the j = 1 contributions in 
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diagrams of type c, e., and f. Again, tlie reason for this is that the combinatoric factors 
are the same. This general behavior has been expUcitly checked for a large number of cases 
and appears to hold for diagrams containing any number of clusters and for clusters of 
any size; in all cases cancelations occur in which the leading behavior after cancelations 
scales as n but not higher. Moreover, it was also explicitly checked to hold for off-diagonal 
matrix elements such as seen in diagram g. of Fig. |l]; these off-diagonal contributions only 
contribute for j = 2 and higher. 

Equation ( |5.5D follows directly from these considerations. Note the underlying ampli- 
tudes associated with clusters such as and T, |_ are computable perturbatively. More- 
over, the number of distinct cluster types which can contribute is fixed by the order of 
perturbation theory /: a cluster of size n requires at least n — 1 gluon exchanges. Thus 
when working at any fixed order I, one needs only amplitudes for clusters of size I + 1 or 
less. It should also be clear why the coefficients in Eq. ( ^.51) for all n > / are universal. 

Next evaluate the derivative in Eq. ( p-5\ ). The right-hand side only depends on r 
through the running of the coupling, which is given by a large Nc RG equation with a /3 
function calculable in perturbation theory: 

— 3 = r > bi{Ncas) (5.6) 

ar ^-^ 

i=l 

where the bi are fixed coefficients. To order {NcttsY'^^, R{n,r) is accurately given by 



i?(n, r) 



I I 



Ci = ^ ^ gmbm'^m+m' ,i (5.7) 
m=l m'=l 
I I 

di = hibmSm+m' ,i 



m=l m' = l 



where gi and hi are the coefficients in Eq. (5.5). Note that Ncas{fJ-'^) is a monotonically 
decreasing function and asymptotes to zero. Thus, setting fi'^ = one can always find 
a value ro such that NcOg evaluated at/x = l/r>l/ro can be made as small as one likes. 
For sufficiently small N^Os the series is dominated by its first term (i = 2) and then clearly 
satisfies Eq. (^). Thus the last condition needed to establish a Hagedorn spectrum is 
satisfied provided standard assumptions about the regime of validity of perturbation theory 
apply. 



6. Discussion 



The argument of the previous four sections demonstrates that a Hagedorn spectrum must 
exist for large Nc QCD provided that renormalization-group improved perturbation theory 
is applicable for correlation functions at sufficiently short distances. 
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One obvious drawback of this approach is its reliance on perturbation theory for the 
correlation functions. Of course, there is a general expectation that perturbation theory 
will be accurate at short distances. Moreover, vacuum condensates yielding power law cor- 
rections to perturbation theory (as one has for example in the QCD sum rule approach to 
phenomenology [p!3[ ) can be incorporated in the present approach in a straightforward way. 
The contributions of the condensate factorize in a manner analogous to the perturbative 
corrections and do not affect the conclusions. Nonetheless, the reliance of perturbation 
theory raises the issue that the perturbative expansion is asymptotic. This appears to 
preclude a straightforward method for strengthening into a rigorous theorem the heuristic 
reasoning used here. However, despite its lack of rigor the approach does give a strong ar- 
gument as to why a Hagedorn spectrum is expected without any reliance on the assumption 
of stringy dynamics. 

The present approach can be generalized easily. With modest changes, it can be used 
to show that mesons and glueballs of any quantum number also have exponentially growing 
spectra. The modification of the argument for mesons with quantum numbers other than 
scalar (and pseudoscalar) involves little more than altering the form of the dispersion 
relation to account for spin. As in the case of spinless mesons, one can insert arbitrary 
numbers of 0± operators between the quark creation and annihilation operators without 
affecting the quantum numbers of the operator. Thus the number of operators grows as 
2" and the argument goes through essentially unchanged. 

There is a small subtlety for the case of glueballs. Consider the case of scalar and 
pseudoscalar glueballs. The natural class of operators to pick are of the form 

Ji,,i,,...,i^ = Ti{Oi,Oi,---OiJ where ,/„ = ± (6.1) 

and the trace is in color space. Because of the cyclic property of the trace these operators 

are not all distinct. For example J+-| is identical to J-| |-. This raises the obvious question 

of how many distinct operators exist for any fixed n. As it happens, the combinatorics are 
nontrivial and there is no simple analytic formula for this for general n. However, it is easy 
to show that when n is a prime number the number of distinct operators is (2" — 2n + 1) /n; 
when n is not prime it is always larger than (2" — 2n-|-l)/n. Thus, the number of operators 
still grows faster than any exponential with base less than 2 and the basic argument goes 
through. 

The argument also applies to other QCD-like theories at large Nc- These include 
pure gauge theory, theories with quarks in non-fundamental representations: the adjoint 
representation, the two-index symmetric representation or the two index antisymmetric 
representation |14|. The latter may be of interest for phenomenological reasons — at A'^^ = 3 



it coincides with QCD. The approach also applies to large Nc gauge theory (with and 
without matter fields) for other gauge groups {eg., SO{Nc))- 

All of the systems with Hagedorn spectra considered above have an unbroken effective 
center symmetry at large Nc- This interesting since an unbroken center symmetry is often 
taken as a signature of confinement. In cases such as pure gauge theory or theories with 
adjoint fermions, the center symmetry is manifest. For the other cases it is an emergent 
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symmetry — while not being a symmetry of the theory at finite Nc, the effects of explicit 
center symmetry breaking becomes small as Nc increases and goes to zero at infinite Nc- 

For QCD with fundamental quarks, center symmetry trivially emerges at large Nc since 
quark loops are suppressed while the pure Yang-Mills sector has an explicit center symme- 
try. The emergence of an effective center symmetry at large Nc depends on nontrivial large 
Nc dynamics in other cases. For example, in theories with quarks in the anti-symmetric 
two-label representation! IC], the emergent center symmetry may be thought of as due to 



an orientifold equivalence at large Nc between this and QCD with adjoint matter. Other 
cases can be seen to have an emergent symmetry due to orbifold considerations [^]. It is 
interesting that although all these cases have an emergent center symmetry, the derivation 
of the Hagedorn spectrum did not rely explicitly on confinement in the sense of an unbro- 
ken center symmetry; it only used confinement in the sense of an absence of non-singlet 
physical states. The question of whether an unbroken center symmetry is required for the 
emergence of a Hagedorn spectrum will be explored in future work. 

The author thanks D.T. Son, L.Ya. Glozman, A. Cherman and M. Shifman for useful 
discussions. He is indebted to R. Pisarski for pointing out the connections to refs. [§, ^, and 
to Michael Cohen for an explanation of the combinatorial aspects of the glueball operators. 
The support of the U.S. Department of Energy is gratefully acknowledged. 

A. The lemma 



Before proving the lemma of Eq. (2.2), it is worth remarking that the logic underlying it 
is quite standard in the lattice QCD community. It was realized long ago that the masses 
of the lowest lying states coupling to a set of currents with fixed quantum numbers are 
extractable as the large r limit of the derivative of the eigenvalues of log ^^f"'^(r)^ . Indeed 
attempts to extract excited hadron masses from numerical lattice studies routinely use 



this approach [18 1 . The lemma follows simply by averaging over all the masses extractable 
this way in principle while working at large Nc (where the mass spectrum is discrete) 
and exploiting the following standard properties of the propagator which hold for positive 
values of r and s. 

N N 9^1og(A(r,s)) ^ 51og(A(r,s)) ^ , 

i)A(r,s)>0 ii) \ . < «m) \. ^ as r ^ oo . A.l 

or^ or 

A formal proof of the lemma starts with an easily derived identity: if ^ {a) is a 
sufficiently smooth matrix valued function of invertible matrices then 



tr {ii{a)-^ii'{a)ii{a)-^ii\a)] - tr {ay^^} " (a)) (A.2) 



tr (log (jl{a) 
d^ 

where the prime indicates differentiation. Suppose that H (a) is Hermitian and positive 
definite for all finite positive a and that ^"(a) is negative definite. Then both traces the 
right-hand side of Eq. ( [A.2| ) can easily be seen to be necessarily positive. This implies that 
for a > , -^tr (logf^(a))) > lim^^oo (-^tr (log (^(a)))). Note that U^{r) is 
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Hermitian. Thus, providing that ^^(r) is positive definite and ^^^{r)J is negative 
definite, 

^ tr (log(tf'5(r))) / ^ tr (log (tf '5(r) 



dr \\S\\ r->oo \ dr 



> hm ^ ^— ^ . (A.3) 




The demonstration that tf^(r) is positive definite and d^^ (tf ^(r)) is negative definite 
follows from the definition of 11 "^(r) in Eq. (|2.l| ) which allows one to rewrite 



r as 



oo 



the sum is over mesons. By construction ^^^^^ is positive definite: it is a positive constant 
times a projection operator on to a single vector. Combining this with property i) of 
Eq. (A.l) implies that ^^{r) is positive definite. Similarly, combining it with property ii) 



implies that (tf^(r)) is negative definite. 
The final step is to show that 



tr (log (tf^lr))) 




which together with Eq. ( [A. 3D establishes the lemma. To do so consider the sum in 
Eq. (A.4). One can construct a new operator by truncating the sum over mesons: 

I?Lnc(0 = E '^^'^^A(r, m|) , (A.6) 

where the mesons included in this sum are in a set Ti with the following properties: i) 
1 1^1 1 — ll'^ll (^-c, there is one state per operator in 5 ); ii) elements of all of the matrices 
for mesons in the set are linearly independent and nonzero; and iii) the mesons 
in the set are the lightest ones possible consistent with previous conditions. Generically 
one expects that in the absence of symmetries this set will simply be the lightest ||5|| 
mesons. Condition iii) of Eq. ( |A.1| ) implies that as r — )• oo , ^"^(r) — )• ^^^^^^.(r) since the 
effects of the states not in T-L make exponentially suppressed contributions which vanish 
as r — )• oo. The exponential suppression of higher mass states at large r also greatly 
simplifies the computation of eigenvectors and eigenvalues of ^^{r) at large r. Recall 
that up to a multiplicative constant the matrices are projections on to a single 

vector Vj. The eigenvectors of ^tounc(^) large r can be constructed from these. Clearly 
vi becomes an eigenvector at large r since the contributions of all mesons in the set except 
the lightest makie a contribution which is exponentially small compared to the leading one. 
The associated eigenvector is A(r, mf)^||f-| jci^^p. One obtains a second eigenvector at 
large r by starting with V2 and projecting out of vi. Higher mass states exponentiate away 
while the contribution from the lightest state is removed by projection; the eigenvalue is 
A(r, m|) X]i!=i 1^2, fcP- One can construct all eigenvectors in a similar manner; the approach 
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is the essentially Gram-Schimdt procedure: the third eigenvector is with V2 and vi 
projected out and so forth. The j*^ eigenvalue is thus A(r, m|) |cj,fcP for all j G Ti. 

Together with condition iii) of Eq.( [A.l| ) this implies that 




(A.7) 



However, Ylj^-H ]\s\\ — Si=i ]\^\ (either H includes the ||5|| lightest states , or due to a 
lack of linear independence or vanishing coefficients, it includes some more massive states). 
Moreover, it has already been shown that tf'5(r) tfg„n^(r). Together with Eq. { \A7i\) 
this implies Eq. ( [A.5| ) and thus establishes the lemma. 
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